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Study of tensor states in the reaction 77 — > tt + 7t tt° 
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Abstract 

A high statistics study of the reaction 77 — ► 7r + 7r~7r° has been performed with the Belle detector 
using a data sample of 26 fb _1 collected at y/s = 10.58 GeV/c 2 . A spin-parity analysis shows 
dominance of the J p = 2 + helicity 2 wave for three-pion invariant masses from 1 to 3 GeV/c 2 . The 
invariant mass distribution exhibits 02(1320), a2(1700) and higher mass enhancements. 

PACS numbers: 12.38.Qk, 12.39.-x, 12.40.Vv, 13.60.Le, 14.40. Cs 
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1. INTRODUCTION 



Three-pion final states of two-photon interactions are restricted to quantum numbers 
suitable for study of resonance formation. The pn channel is known to be dominated by 
the formation of a 2 (1320) in the helicity 2 state [1-10]. The a2(1320) is a ground state of 
isospin 1 3 P 2 QQ mesons. Observations of higher mass resonances have been reported [7-10]. 
Study of higher mass states is important for the assignment of nonet members and for the 
understanding of confinement in the quark model [11-15]. The L3 collaboration reported the 
a2(1700) in pn and f%n decay modes mainly in helicity 2 states Tol . The Belle observation 
of a resonance in 77 — > K + K~ [16] at 1737 MeV/c 2 is listed in the Particle Data Group 
(PDG) [ItJ and attributed to the a 2 (1700). Experiments on np collisions [18-20] and jp 
photo-production 21, 2^1 have also observed a three-pion resonance near 1.8 GeV/c 2 . The 
Crystal Barrel collaboration reported the ci2(1700) in the n°rj mode in pp collisions j2j| 24] . 
The 02(2100) and a 2 (2280) in the /2 71 " mode were reported in |25j . 



We present a sp in-p arity analysis of three-pion events using data collected with the Belle 
detector at KEKB [26j . Resonance formation is investigated in the pir and /2 71 " modes, includ- 
ing interference between them. The data sample was taken with similar trigger conditions 
during 2000 and 2001 at a center-of-mass energy y/s = 10.58 GeV. The corresponding total 
integrated luminosity is 26.0 fb _1 . We first describe the theoretical formulae for spin-parity 
analysis of three-pion events. The spin dependence is investigated using the distribution 
of the vector normal to the three-pion decay plane. Upper limits on two-photon radiative 
widths are determined for 7r (1300), 7r 2 (1670) and a 4 (2040). A neural network method is 
applied to enhance tensor state selection. The helicity states of tensor resonances are inves- 
tigated using cos 9 distributions of final state pions. Background contamination is evaluated 
by comparing p 2 (37r) distributions for data and Monte Carlo generated two-photon interac- 
tions. The three-pion and di-pion mass spectra are investigated. The two-photon radiative 
width of the a 2 (1320) is measured and compared to the PDG world average. The mass region 
above 1.5 GeV/c 2 is examined for the a2(1700) and new states. The two-photon radiative 
widths are determined. Interference between tensor states and pn and /27r decay modes is 
investigated and the coupling amplitudes and phase angles are also determined. 



2. THEORETICAL FORMULAE AND MONTE CARLO 



The three-pion final state in two-photon interactions is expressed as the composition of 
production of di-pion isobars, In, and the decay of / —>■ nn. Possible di-pion isobars are fo, 
p ± and / 2 . The assignment of spin-parity {J p ) and orbital angular momentum (L) are listed 
in Table HI Gauge invariance and Bose symmetry forbid J p = 1^ and 3~. Only the helicity 
A = state is allowed for J p = 0" and 2", and A = 2 for J p = 3+ [27]. The helicitv 
fraction of a tensor state is predicted to be zero by the non-relativistic quark model 28[ 2t| . 
The static quark model also predicts a small value (1/7) (3oT |. 

The cross section for resonance formation in 77 — > In is given by 



da 77 = 2vr(2 J + l)f 77 ^ R d 
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TABLE I: Spin-parity of the three-pion final state in two-photon interactions. Listed are the 
possible di-pion isobars (/q, p or / 2 ), helicity (A) states allowed, and the lowest orbital angular 
momentum (L) between the di-pion isobar and the third pion. 

where T 77 is the two-photon radiative width of the ground state resonance of spin J and 
Rj z is the probability of helicity J z . BWj = l/(s — mf + im^i) is the Breit-Wigner term for 
two-photon invariant mass (W^ = \fs) for a resonance of mass and width 1~V The decay 
modes are denoted by Dj z (I). Interference with the ground state is given by the coupling 
amplitude a« and phase angle fa, where the radiative width is afT 77 . The decay amplitude 
is given by 

= BW(p + )T J *(p+,7r) +BMj>-)T J '(j>-,n) + &e^BW(/ 2 )T J *(/ 2 , vr) (2) 

i 

for the Breit-Wigner terms of p and / 2 isobars with interference expressed by the amplitude 
£i and phase angle ipi- The spin dependence is given by 



T J *(/,7T) = 32VT 2 ( mR T R m 27T T 2n ^^) Y. CJ d-rn,l, m YL Z ~ m (^A2«)Yr(6.An), (3) 
V P2nPir J m 

where the subscript 2tt denotes parameters of the the di-pion (p ± or / 2 ) isobar with spin / and 
third component m. The subscript tc denotes a pion from the di-pion decay. The spherical 
harmonics are multiplied by the Clebsch-Gordan coefficients C J ^\_ mlm . They describe the 
angular distribution of the di-pions (with invariant mass m 27r = v/s^, momentum p 2?r , polar 
angle # 27r , and azimuthal angle 2yr ) in the resonance rest frame and the pions from the 
di-pion decay (with p n , 6 n , </>„.) in the di-pion rest frame. The 6 and </> angles in both frames 
are determined with respect to the incident 77 direction which is approximated by the e + e~ 
beam direction in the center-of-mass frame. 

The Monte Carlo is prepared for two-photon interactions in the asymmetric beam con- 
figuration (8 GeV e~ on 3.5 GeV e + ) at Belle. The incident photon flux is calculated with 
the two-photon luminosity function in j3jj. The colliding photons are generated with the 
p-pole form factor of the vector meson dominance model (VDM). The photon flux has been 
compared to the QED calculation of the DIAG36 program js^]. The systematic uncertainty 
on the photon flux simulation is about 2%. 

The detector response is simulated with GEANT3 [33|. The event trigger is simulated 
using trigger inputs from the Belle tracking system and electromagenetic calorimeter. 

3. EVENT SELECTION 

Events are pre-selected for two-photon interactions with the energy sum of final state 
particles below 5 GeV and the scattered electrons not observed. Charged pions are detected 
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by the central drift chamber (CDC) and the silicon vertex detector (SVD) in a 1.5 T magnetic 
field. The coverage in polar angle extends from 17° to 150° in the laboratory frame. Outside 
the CDC are the Aerogel Cerenkov Counter system (ACC) and the Time Of Flight (TOF) 
system. The photons from 7r° decays are detected by the electromagnetic calorimeter (ECL) 
comprised of CsI(Tl) crystals covering the same angular region as the CDC. A detailed 
description of the detectors and their performance is given in reference |34| . 

The three-pion event trigger is composed of a CDC two-track trigger, an ECL energy 
trigger, and combinations with looser constraints on the track opening angles and energy 
thresholds. In event reconstruction, charged pion tracks are required to have transverse 
momentum larger than 300 MeV/ c. The purity of the charged pions is enhanced by rejecting 
electrons using a likelihood function based on the ratio of ECL energy to track momentum 
and dE/dx of the track in the CDC |35j. Charged kaons are rejected by cuts on a probability 
function based on ACC, TOF and CDC measurements. ECL and beam background are 
suppressed by requiring a maximum distance to the interaction point (IP) of 2.5 mm in the 
r-(f) plane and 50 mm along z (the beam direction). A photon is identified as an isolated 
ECL cluster with no matching track. Partially contained photons in the ECL and beam 
background are eliminated by requiring an energy threshold E 1 above 100 MeV. 

Three-pion candidates are selected by requiring two oppositely charged pions and a pho- 
ton pair with invariant mass within three sigma of the 7r° mass. Background events are 
suppressed by a tight cut on the three-pion transverse momentum and a high energy thresh- 
old for photons from 7r° decay. Additional selection criteria are chosen for the event topol- 
ogy. Distributions of charged pions and photons are simulated for the partial waves and 
decay modes to be investigated. In the three-pion rest frame, the photons from 7r° de- 
cay are expected to have a large opening angle (^(77)) and to make large angles with the 
charged pions (^2(777)). Background is expected from inclusive events or low multiplicity 
events contaminated with random photons. These events populate the low opening angle 
region and can be eliminated by requiring min(f2(7r7) + ^(77)) > 1 radian. Background 
events with a di-pion in the low mass region outside the Dalitz contour are excluded by 
m(7r + 7r°) +m(7r~7r°) > 1 GeV/c 2 . Contamination by 77 — > 7t + tt~ events in coincidence with 
a fake tc° is suppressed by requiring p 2 (7r + 7r~) > 0.01 GeV 2 /c 2 , as the transverse momentum 
of the 7i + tt~ pair is balanced in direct production. 



4. SPIN-PARITY ANALYSIS 

The spin-parity of a three-pion event is characterized by the spherical harmonics in Eq. |21 
The dependence is explicitly seen in the angular distribution of the final state pions and of 
the vector normal to the three-pion decay plane. The spin-parity of the data sample is 
identified using the A parameter, which is the squared magnitude of the normal vector, 
scaled by the maximum available kinetic energy [J 0] , 



A 



N/Q 



(4) 



where N = p 7T + x p n - is evaluated in the three-pion rest frame and Q is the rest energy 
difference between the resonance and the three pions. The normal vector N covers a wider 
cos# range than the detector acceptance, providing better discrimination than the the cos 6 
values of individual pions. 
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Goodness-of-fit tests to the A distributions in Fig. ^ were performed in four invariant mass 
intervals. The event sample is selected with p 2 (37r) < 0.001 GeV 2 /c 2 and E 1 > 150 MeV. 
The mass intervals were chosen to select the 02 (1320), the 02 (1700) and possible higher 
mass states. Monte Carlo distributions were generated for resonances of each partial wave 
hypothesis investigated. The generated resonance masses and widths correspond to the 
enhancements observed in data (to be discussed in Section 7). Interference between the pn 
and f27r° channels is included using £ = 0.91 and ip = 150°. 

The a2(1320) is dominant in the mass range 1.0 GeV/c 2 < m(37r) < 1.5 GeV/c 2 (Fig.[TJi). 
The A distribution is consistent with the J p = 2 + expectation in the pir decay mode. The 
test results for higher mass regions are also consistent with the J p = 2 + expectations for 
the and fiit modes. The \ 2 values obtained are listed in Table |H] The good agreement 
seen in the low A region indicates little contamination by phase space background. The 
helicity state is not distinguished by the A distribution, but is better determined by the 
cos# distribution of final state pions (to be discussed in Section 6). 

Contributions of known non-tensor resonances, 7r (1300)(l~), 7r 2 (1670)(2~), and 
a4(2040) (4 + ) are estimated by binned maximum likelihood fits to the A distributions. The 
Monte Carlo distributions are generated for a test resonance of spin-parity J p and a tensor 
resonance (of helicity-2) of mass and width listed in Table HP for the chosen mass range. The 
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FIG. 1: A distributions in comparison with Monte Carlo of various spin-parity hypotheses. Monte 
Carlo distributions in each mass interval are simulated for resonances of mass and width listed in 
Table HH 
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TABLE II: x 2 / n df of goodness-of-fit tests to the A distributions of spin-parity and helicity states. 
Monte Carlo simulations were performed with the resonance masses and widths listed. 

sum of the distributions, with the bin contents 

N MC = f . N MC {jP) + (1 _ /)jV MC (2+)j (5) 

is fitted for the fraction /. 

The contribution of the 7r (1300) is found to be negligible in the A distribution of m(37r) < 
1.50 GeV/c 2 for the pir and it + it~it decay modes. The result / = 0.0±1.1±1.0% is obtained 
for both decay modes. The 0~ wave is easily distinguishable from the tensor wave and there- 
fore a small systematic uncertainty (1.0%) is estimated for the dependence on mass resolution 
and selection efficiency. The expectation for the 7r (1300) is simulated with the resonance 
parameters from the PDG, m = 1300 ± 100 MeV/c 2 and width in the range T = 200 — 600 
MeV/c 2 . The number of events estimated for r 77 (7r (1300))i3(7r (1300) -> pvr) = 1 eV is 
n mc = 3 2 4 ± 0.32 for the pn mode, and 2.93 ± 0.29 for the tt+tt-tt mode. The errors 
are estimated from the uncertainties on mass and width. The upper limit is derived from 
the 7To(1300) event fraction found in the fit and the quadratic sum of statistical and sys- 
tematic errors (1.49%). It corresponds to 173 of the total 11654 data events. Applying 
Poisson statistics, the upper limit at the 90% confidence level (CL) is 191 events, corre- 
sponding to two-photon radiative widths r 77 (7r (1300))£>(7ro(1300) — > pn) < 65 eV and 
r 77 (7r (1300))£(7r (1300) -> tt+tt-tt ) < 72 eV, respectively. The N MC used in the calcula- 
tion is reduced by one standard deviation to account for the uncertainties of the 7r (1300) 
mass and width. 

Likewise, the contribution of 7r 2 (1670) decaying into p-n and / 2 7r is determined with the 
A distribution of 8836 data events satisfying 1.5 GeV/c 2 < m(37r) < 1.85 GeV/c 2 . The 
event fraction obtained for the 7r 2 (1670) is f — 0.0 ± 1.2 ± 2.0%. The expected number 
of events for r 77 (vr 2 (1670))i3(7r 2 (1670) -> pn, / 2 tt) = 1 eV is 19.3 ± 0.4 events for tt 2 (1670) 
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with m = 1672.3 ± 3.2 MeV/c 2 and T = 259 ± 9 MeV/c 2 . The upper limit obtained is 
r 77 (7r 2 (1670))i3(7r 2 (1670) -> pn, f 2 n) < 12 eV (90% CL). 

The contribution of a4(2040) is determined for helicity and 2 states with the A distribu- 
tion of 4400 data events with 1.85 GeV/c 2 < m(37r) < 2.1 GeV/c 2 . The fractions obtained 
in this case are also negligible, with / = 0.0 ± 4.5 ± 4.0% for both helicity states. The 
numbers of events expected for r 77 (a 4 (2040))£(a 4 (2040) -> pn, / 2 tt) = 1 eV is 48.2 ± 2.2 
and 46.4 ± 2.1 for helicity and 2, respectively. The upper limit on the helicity 2 wave is 
found to be r 77 (a 4 (2040))£(a 4 (2040) -»■ pn, / 2 tt) < 6.5 eV (90% CL). 



5. SELECTION OF TENSOR STATES 



5.1. Neural network method 



A feed-forward neural network is employed to improve the selection purity of the chosen 
partial waves and decay modes. The neural network contains seven input nodes and multiple 
outputs, each corresponding to a partial wave. The input variables include the di-pion 
masses and the cos# values of the final state pions and of the vector normal to the decay 
plane. The three-pion events are best described by a tensor partial wave in the pn and f 2 n 
modes. The neural network is therefore trained to discriminate helicity 0, helicity 2, and 
background (approximated by phase space). The product of the outputs of three nodes Oj, 
for J p (\) = 2 + (2), 2 + (0), and phase space, gives the neural network weight 

m = o i -(i-o j )-(i-o h ), (6) 

where % is the wave chosen to be identified against the other two. Distributions of NN for 
helicity 2 are shown in Fig. El The data sample used has p 2 (37r) < 0.0005 GeV 2 /c 2 and 
photon energy E 1 > 180 MeV. The sample is divided into masses below and above 1.5 
GeV/c 2 for comparison of spectra with events dominated by the a 2 (1320) and those in the 
higher mass region. The NN spectra are in good agreement with the helicity 2 (hatched area) 
predictions. Note that the detector acceptance is higher for the helicity sample (dashed 
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FIG. 2: Neural network output for J p = 2 + helicity 2 wave. The Monte Carlo samples are 
normalized to fixed radiative widths for comparison of selection efficiencies of different partial 
waves. r 77 £> = 0.7 keV is applied (for a 2 (1320)) in (a), and r 77 £> = 0.5 keV in (b), respectively. 
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line). The selection efficiency for the phase space sample is comparable with that for the 
tensor states; however, the event rate differs by the spin factor (2J+ 1) in the cross section. 
With a cut applied at NN > 0.2, phase space events at lower mass are suppressed by 80%, 
leaving about equal numbers of events in the helicity and 2 states. The a<i (1320) selection 
efficiencies become 0.023%, 0.053% and 0.054% for phase space, tensor helicity and tensor 
helicity 2 decays, respectively. 

The event selection cuts were optimized for selection purity. Because events from three- 
pion resonances are densely distributed at pt(37r) near zero, a p\ (3ir) cut is effective in 
improving selection purity. A relatively high photon energy threshold suppresses beam 
related background and an NN threshold suppresses random three-pion background. A 
data sample enriched in events from tensor states is selected with p\ (3ir) < 0.0005 GeV 2 /c 2 , 
£ 7 > 180 MeV, and NN(2 + (2)) > 0.2. The selection efficiency is evaluated with Monte Carlo 
events. Figure E3 shows the generated and selected mass spectra for the J p = 2 + helicity 2 
wave in the pn and /2 71 " modes for three resonances at 1320, 1750 and 1950 MeV/c 2 . The 
selection efficiency as a function of m(37r) is a smooth curve. The shape does not change 
much as the cut values are varied, which helps in reducing the systematic uncertainty. 
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FIG. 3: Estimated selection efficiency for the J p = 2 + helicity 2 partial wave. Distributions are 
shown for Monte Carlo generation and fully simulated and selected events. The efficiency, shown 
as points with error bars, is fitted to a polynomial function. 



5.2. Background template 

Background contamination is investigated with rejected events. Their p 2 (37r) distribution 
is approximately linear, which differs very much from two-photon events with a p-pole form 
factor peaking at zero. The background fraction is estimated from a fit to the p 2 (37r) 
distribution of Monte Carlo two-photon events combined with a linear background function. 

The sample is divided into 100 MeV/c 2 mass intervals in m(37r); some of the fits to pi 
distributions are illustrated in Fig. |3J The background fractions are then obtained for a 
chosen p 2 (37r) threshold, and are parameterized as a function of m(37r). The background 
distribution thus obtained is a smooth curve corresponding to 15% of the final sample of 
9291 events (Section 5.1) in the three-pion mass range up to 3 GeV/c 2 . 



10 




J 0.004 0.004 0.004 0.004 

p2(37t) GeV 2 /c 2 

FIG. 4: Some pi (3n) distributions in 100 MeV/c 2 mass intervals. Each distribution is fitted to the 
Monte Carlo simulation and a linear background. 
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FIG. 5: cos distributions of final state pions. Solid lines are the fits of helicity fractions. 



6. HELICITY AMPLITUDE 

The helicity fractions of the selected three-pion events are investigated by examining the 
polar angle distributions of pions. The cos 9 distributions for the of 7r° and the difference in 
cos 9 of 7r + and tt~ are shown in Fig. in four invariant mass intervals, corresponding to 
02(1320) and possible higher mass states. The two helicity states are clearly distinguishable 
in cos#(7r°) with helicity 2 (0) events distributed in the forward (central) region, and vice 
versa for the distributions of cos9(tt + ) — cos#(7r~). The asymmetric shape observed in 
cos6*(7r°) distributions is due to the Belle configuration, which has asymmetric acceptance 
along and opposite to the boost direction. 
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?7i(37r) mass range 



background (fixed) % Helicity-2 % x 2 /ndf 



1.0 - 1.5 GeV/c 2 
1.5 - 1.85 GeV/c 2 
1.85 - 2.1 GeV/c 2 

2.1 - 2.4 GeV/c 2 



0.8 100 ±2 ±5 40/35 

10 95 ± 2 ± 5 96/37 

13 100±2±5 131/36 

30 82 ± 2 ± 5 39/37 



TABLE III: Helicity 2 fractions of the fits to cos 9 distributions. Background fractions estimated 
in Section 5.2 are fixed in the fits. 

The helicity fractions are obtained by fits to the cos 9 distributions for the two helicity 
states and background. The fits, shown as histograms in Fig. are consistent with the 
hypothesis of helicity 2 dominance in all the mass intervals. The fits with the background 
fractions fixed to the estimates in Section 5.2 are listed in Table InTl Background is assumed 
to be randomly distributed in three-pion phase space. Statistically compatible results are 
obtained with background fractions set to zero. In the low mass region, where a 2 (1320) is 
dominant, the cos# distributions are consistent with pure helicity 2. In the higher mass 
region the helicity state is also predominately helicity 2. The Monte Carlo includes multiple 
tensor states in the pit and ftft channels. Note the discrepancy in cos 0(ir°) for 1.85 GeV/c 2 < 
m(37r) < 2.1 GeV/c 2 , where the Monte Carlo suggests more energetic final state pions and 
therefore a more symmetric distribution. Systematic uncertainties are estimated by varying 
m(3ir) intervals, selection cuts and background levels. The deviation in helicity fractions is 
evaluated and the systematic errors estimated are listed in Table ITTT1 

7. MASS SPECTRUM OF TENSOR STATES 
7.1. Resonance parameters 

The invariant mass spectrum of three-pion events is shown in Fig. El The a2(1320) is 
dominant in the low mass region where the mass spectrum is truncated by trigger thresholds 
and selection cuts. The a2(1320) radiative width is determined by the number of events 
observed in the mass region below 1.5 GeV/c 2 . The Monte Carlo prediction has included 
full detector simulation with the 09.(1320) resonance parameters fixed to the PDG values 
(m = 1319 MeV/c 2 ,T = 105 MeV/c 2 ) 0. The background estimated in Section 5.2 is 
included. By scaling the Monte Carlo expectation for a2(1320) shown by the solid line in 
Fig. IHl the two photon radiative width obtained is 



It is in good agreement with the PDG world average. The systematic uncertainty is esti- 
mated to be 11%. It is attributed to the trigger efficiency and event selection (7%) and the 
background and contribution from higher mass states (8%). 

The spectrum in the higher mass re gion shows a structure too broad to be due only to the 
o 2 (1700), with width about 200 MeV The quick rise at 1.5 GeV/c 2 and the dip at 1.8 
GeV/c 2 are indications of overlapping resonances. The interference of tensor resonances is 
parameterized in Eq. ^ by the coupling amplitudes and phase angles relative to the ground 
state a 2 (1320). The fit for the resonance parameters was conducted with the mass spectrum 



r 77 (a 2 (1320)) = 0.99 ± 0.03 ± 0.11 keV. 
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Data 

a 2 (1320) MC 
Solution-1 
Solution-2 
Bkgd 




m(3jt), GeV/c' 

FIG. 6: Three-pion mass spectrum. The fit to a2(1320) in the mass region below 1.5 GeV/c 2 is 
shown by the solid line. The two solutions of the fit to interference of tensor states are also shown. 
The dashed-dotted line is the estimated background. 



generated by Monte Carlo simulation. The ai (1320), serving as the reference ground state, 
is generated with the resonance parameters given by the PDG and the resulting spectrum 
is normalized to data in the mass region below 1.5 GeV/c 2 . 

The fit giving satisfactory agreement with the data is achieved with four resonances, which 
results in eight sets of solutions. The interference among the tensor states has the effect of 
shifting resonance mass positions and shapes. The coupling to the phase angles is periodic 
with period 2tc. In each half period, one solution is obtained with the amplitudes of higher 
mass states smaller than 1 and closely grouped phase angles. These solutions are shown by 
the dashed and dotted lines in Fig.Elfor Solution-1 and 2, respectively, with x 2 /ndf = 45/50 
and 46/50 in the mass range up to 2.4 GeV/c 2 . The parameters obtained are listed in 
Table IIVI The solutions with interference amplitudes larger than 1 (Solution^ > 1), 
Table IIV|) correspond to destructive interference with radially excited states having larger 
cross sections than the ground state (a2(1320)), which is unlikely for light-quark mesons. 

Correlations between the fitted phase angles and resonance parameters are significant. 
The resonance mass positions are shifted to values higher than in the fit without interference. 
The largest correlation coefficient is between (pi and m(1700) for a2(1700). The values 
obtained are -0.72 (-0.88) in Solution-1 (2), respectively. The width of the resonance is 
positively correlated with the interference amplitude. The relative amplitude «2 has a large 
correlation coefficient with the width of the resonance at 1950 MeV/c 2 . The values are 0.74 
(0.75) in Solution-1 (2), respectively. The interference of overlapping resonances leads to 
strong correlations and to shifts in two-photon radiative widths. The radiative width of 
a2(1700) in three-pion mode is found using r 77 (1700)£>(37r) = a\ ■ r ry (1320)£>(p7r) to be 
0.10 ± 0.02 ± 0.02 keV. 

Systematic uncertainties are evaluated including the effects of interference on selection 
efficiency and mass resolution. The mass resolution is smeared in the reconstruction of 
7T° and charged pion tracks. The smearing effect on mass is estimated by Monte Carlo 
and is found to be approximately equivalent to convolution with a Gaussian of width less 
than 20 MeV/c 2 . The Monte Carlo spectrum is slightly shifted with respect to data and the 
magnitude of the shift is examined as a function of the trigger and event selection thresholds. 
The precision of the resonance mass measurement is estimated to be better than 8 MeV/c 2 . 
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The uncertainty in the resonance width is due to the smearing of the invariant mass 
measurement. It is estimated to be 10 MeV/c 2 by comparing the Pi(37r) of data and Monte 
Carlo and by comparing the simulated and measured widths of the 02 (1320). The systematic 
error from interference amplitudes is 9%. It is estimated using the uncertainties of selection 
efficiencies and backgrounds. The interference phase angles are strongly correlated with the 
resonances masses. The systematic error on masses corresponds to an error of 8 degrees 
for the phase angles in Solution- 1 and 2. In addition, the phase angles are affected by the 
decay modes, and the errors are evaluated in the fits to be 8 degrees. The combined error 
is estimated to be 12 degrees. The fits having large amplitudes (ctj > 1) are more sensitive 
to phase angle (smaller statistical error) and are less dependent on the errors in the masses. 
The systematic uncertainty of the phase angles is thus smaller. 



Solution- 1 


mass MeV/c 2 width MeV/c 2 amplitude (a's) phase (</>'s) deg. 


X 2 /ndf 


a 2 (1700) 
BW(1950) 
BW(2U0) 


1769 ±10 ±8 270±7± 10 0.37 ± 0.02 ± 0.03 154±6±12 
1948 ±4± 8 291±6±10 0.61 ± 0.02 ± 0.06 143 ± 5 ± 12 
2146 ± 12 ± 8 358 ± 26 ± 10 0.40 ± 0.03 ± 0.04 139 ± 4 ± 12 


45/50 


Solution-2 


mass MeV/c 2 width MeV/c 2 amplitude (a's) phase (</>'s) deg. 




a 2 (1700) 
BW(1950) 
BW(2U0) 


1758 ± 13 ± 8 269 ± 10 ± 10 0.37 ± 0.03 ± 0.03 221 ± 7 ± 12 
1949 ± 6 ± 8 324 ± 14 ± 10 0.71 ± 0.03 ± 0.06 220 ± 6 ± 12 
2161 ± 17 ± 8 342 ± 22 ± 10 0.44 ± 0.03 ± 0.04 221 ± 6 ± 12 


46/50 




Solution^, > 1) 


mass MeV/c 2 width MeV/c 2 amplitude (a's) phase (4>'s) deg. 


X 2 /ndf 


a 2 (1700) 
J3W(1950) 
BW(2U0) 


1768±5±8 273±8±10 0.50 ± 0.02 ± 0.05 126 ± 2 ± 8 
1949±3±8 280±4±10 1.54 ± 0.02 ± 0.14 99 ± 1 ± 8 
2138 ±4± 8 382±5±10 2.03 ± 0.03 ± 0.18 236 ± 1 ± 8 


48/50 


a 2 (1700) 
.BW(1950) 
BW(2U0) 


1757 ± 6 ± 8 279 ± 12 ± 10 0.54 ± 0.02 ± 0.05 207 ± 3 ± 8 
1954±3±8 289±5±10 1.64 ± 0.02 ± 0.15 163 ± 1 ± 8 
2136 ±4± 8 374±4±10 2.23 ± 0.03 ± 0.20 313 ± 1 ± 8 


49/50 


a 2 (1700) 
BW (1950) 
BW(2U0) 


1750 ±4± 8 271±3±10 1.45 ± 0.02 ± 0.13 102 ± 1 ± 8 
1960±4±8 267±5±10 1.02 ± 0.02 ± 0.09 237 ± 1 ± 8 
2138 ± 9 ± 8 375 ± 14 ± 10 0.64 ± 0.02 ± 0.06 212 ± 2 ± 8 


49/50 


a 2 (1700) 
SW(1950) 
BW(2U0) 


1759 ±4± 8 269±4±10 1.70 ± 0.02 ± 0.15 161 ± 1 ± 8 
1963 ±4± 8 263±5±10 1.31 ± 0.02 ± 0.12 305 ± 2 ± 8 
2136 ± 9 ± 8 360 ± 12 ± 10 0.76 ± 0.03 ± 0.07 290 ± 3 ± 8 


51/50 


a 2 (1700) 
BW(1950) 
BW (2140) 


1745±2±8 268±2±10 1.88 ± 0.02 ± 0.17 81 ± 1 ± 8 
1947±2±8 275±2± 10 2.78 ± 0.02 ± 0.25 188 ± 1 ± 8 
2143±2±8 372±4± 10 2.29 ± 0.02 ± 0.21 305 ± 1 ± 8 


49/50 


a 2 (1700) 
BW (1950) 
BW(2U0) 


1742 ±2±8 263±2±10 1.94 ± 0.02 ± 0.17 143 ± 1 ± 8 
1944±2±8 263±2±10 2.81 ± 0.02 ± 0.25 257 ± 1 ± 8 
2148±2±8 378±5± 10 2.42 ± 0.03 ± 0.22 20 ± 1 ± 8 


51/50 



TABLE IV: Resonance parameters of the fit to the invariant mass spectrum for tensor states with 
interference expressed in Eq. ^ 
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7.2. Intermediate decay channels 



The presence of pit and /27r intermediate states is demonstrated by the di-pion invariant 
mass spectra shown in Fig. [7| The a2(1320) decays into pn only. The m(7r ± 7r°) spectrum 
for m(3ir) < 1.5 GeV/c 2 (Fig. E^) is consistent with simulated pir decay and the m(7r + 7r~) 
spectrum from phase space. The f% is present in the m(7r + 7r~) distributions of m(37r) > 1.5 
GeV/c 2 . The interference between pit and modes is formulated with an amplitude and 
a phase angle in Eq. The tensor resonance parameters have little effect on the di-pion 
invariant mass spectra and the final state pion distributions. Therefore the decay parameters 
are determined independently. 

The interference term has the effect of pulling di-pion mass peaks for p and f^, near their 
nominal mass positions for ip = 180°, to give destructive interference and flatter spectra 
as ip approaches zero. The interference parameters are determined by a least-squares fit 
to the two di-pion mass spectra. The Monte Carlo expectation includes the resonances in 
Section 7.2. To account for the correlation between resonances, the fit was performed for 
the combined sample of 02(1700) and the resonance at 1950 MeV/c 2 with 1.5 GeV/c 2 < 
m(37r) < 2.1 GeV/c 2 . The background included is a phase space distribution scaled to the 
fraction estimated in Section 5.2. The fit obtained with x 2 /ndf = 92/51 is illustrated in 
Fig. [TJd. The parameters obtained are listed in Table The interference amplitudes and 
phase angles are similar for the two resonances. The fit was also conducted for the spectra in 
narrower m(37r) intervals corresponding to the dominance of a<i (1700) and the resonance at 
1950 MeV/c 2 . The Monte Carlo was simulated with the event fractions of resonances varied. 
The results are statistically compatible and are used to estimate systematic uncertainty. 

The event fraction attributed to the resonance at 2140 MeV/c 2 is only about a quarter 






mCTiV) GeV m(7tV) GeV 



FIG. 7: Di-pion mass spectra and the fits to pir and /27r modes. The interference parameters are 
listed in Table IVl 
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amplitude (£'s) phase (V>'s) deg. 



a 2 (1700) 
BW(1950) 
BW(2U0) 



0.92 ±0.10 ±0.08 151±4±12 



0.91 ±0.10 ±0.12 149±4±20 



1.0 ±0.20 ±0.30 145 ±10 ±30 



TABLE V: Interference parameters for pn and /2 71 " decay modes of the fit to the di-pion mass 
spectra. 

of the total events in the mass range 2.1 GeV/c 2 < m(3ir) < 2.4 GeV/c 2 . The interference 
parameters are again determined by fitting the di-pion mass spectra. The x 2 / n df = 63/59 
is obtained in the fit (Fig. Et). The parameters obtained are also listed in Table IVl 

The large overlap of events in adjacent tensor states leads to large systematic uncertainty. 
To estimate the resulting uncertainties, fits to the di-pion mass spectra of events selected in 
wider mass intervals of m(37c) were carried out. For the resonance at 2140 MeV/c 2 , the error 
is up to 30% for the amplitude and 20% for the angle. The estimated errors are included in 
Table El 

8. CONCLUSION 

In the region of low three pion mass the reaction 77 — ► 7r + 7r~7r° is dominated by the 
a2(1320) in the helicity 2 state. Events in higher mass regions are also found to be domi- 
nated by the spin-parity J p = 2 + helicity 2 partial wave. The contribution of 7ro(1300) is 
negligible. The upper limit is found to be r 77 (7r (1300))£(7r (1300) -»■ tt+tt-tt ) < 72 eV 
(90% C.L). The contributions of the 7r 2 (1670) and a 4 (2040) resonances are also negligible. 
The upper limits are determined to be r 77 (7r 2 (1670))£>(7r 2 (1670) — ► pn, f 2 n) < 12 eV and 
r 77 (a 4 (2040))£(a 4 (2040) -> pvr, / 2 tt) < 6.5 eV (90% CL), respectively. 

The observed radiative width of the a 2 (1320) is consistent with the PDG world aver- 
age. The mass spectrum in the region above the a 2 (1320) consists of a broad distribution 
attributed to the a 2 (1700) and is best interpreted with two more radially excited tensor 
states at 1950 MeV/c 2 and 2140 MeV/c 2 decaying to pir and / 2 7r. The mass spectrum is well 
represented by interfering tensor states with coupling amplitudes and phase angles relative 
to a 2 (1320). The amplitude obtained for a 2 (1700) is 0.37, corresponding to the radiative 
width of r 77 (a 2 (1700))£(a 2 (1700) -> pvr, / 2 tt) = 0.10 ± 0.02 ± 0.02 keV. 

The higher mass tensor states observed here may be compared to the theoretical predic- 
tions of the relativistic quark model [ill EH, • The radiative widths obtained in this study 
are compatible with the relativistic quark model calculations. 
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